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SUMMARY 

The performance of three acoustic boundary condition formulations is investigated. The effect of implementation 
differences is also studied. Details of all implementations are given. Results are shown for the acoustic field of a monopole 
in a uniform freestream. 


INTRODUCTION 

Accurate prediction of flow fluctuations representing sound is the goal of Computational Aeroacoustics (CAA) (see, 
for instance, Ref. 1). To do this, a grid is constructed which covers the region of interest, and the discretized governing 
equations are solved on this grid. The computational domain is usually finite, and boundary conditions must be imposed at 
the edges of the grid. These boundary conditions can generate spurious fluctuations that render the computed solution 
entirely unacceptable. 

Most of the classical boundary treatments used in computational fluid dynamics are concerned only with the steady- 
flow solution. Several suggestions for boundary conditions for unsteady flows have recently been proposed. These 
proposals can be classified in three categories: (I) quasi one-dimensional characteristics, (II) decomposition of the solution 
into Fourier modes, and (ID) asymptotic analysis of the governing equations for large distances. 

The present work is concerned with evaluating these types of boundary conditions by comparing the numerical results 
obtained with various boundary conditions to a known analytical solution. The effect of implementation of each boundary 
condition on its performance is also investigated. 

GOVERNING EQUATIONS 

The equations to be solved are the Linearized Euler Equations, which may be written in cylindrical coordinates as: 

M + ^ + + = s (1) 

3t 3x r 3r r 30 
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where: 


and 


Q = 


p 

U 

V 

w 

e 


( 2 ) 


F = 


p + 2uU — pU 2 
uV + vU-pUV 
uW + wU-pUW 
(p + e)U + (u - pU)E 


(3) 


G = <1 


uV + vU-pUV 
p + 2vV — pV 2 
vW + wV-pVW 
(p + e)V + (v-pV)E 


(4) 


H = 


uW + wU - pUW 
vW + wV-pVW 

p + 2wW-pW 2 
(p + e)W + (w-pW)E 


(5) 


S = 


0 

0 

p + 2wW-pW 2 
vW + wV-pVW 
0 


( 6 ) 


Here 


' = (y-l)>[e- 


uU + vV + wW|+^p|u 


1 x/tt2 +v 2 +w 2 


(7) 


and 


P = (y-1)p 


E-^|u 2 +V 2 +W 2 j 


( 8 ) 


(p. u, v, w, e) = [p’ , (pu)' , (pv)' , (pw)' , (pe)’ ]. 

In this notation, U is the axial mean velocity, V is the radial mean velocity, and W is the azimuthal mean velocity. 


( 9 ) 


2 


NUMERICAL SCHEME 


The code is a modified split MacCormack solver, which is second order accurate in time and fourth order accurate in 
space. This extension of the MacCormack scheme is known as the 2-4 scheme, and was developed by Gottleib and 
Turkel. 2 This scheme has been used successfully before by Farouk et. al. 3 and Ragab and Sheen 4 for studying nonlinear 
instability problems in plane shear layers. Sankar, Reddy, and Hariharan 5 have evaluated this scheme for aeroacoustics 
applications. The solution procedure is as follows: 

Equation (1) may be rewritten in operator form as: 


q n + 1 =L rx q n (10) 

In words, given the flow variables q at a time level ‘n\ the code computes derivatives in both the radial and axial 
directions and uses these quantities to advance the flow field one time step, to the ‘n+1’ time level. 

In the present code, the operator is split into separate radial and axial contributions: 

q n + 1 =L r L x q n (11) 


Each operator solves a portion of Eq. (1): 




3Q , 1 9(rG) 

3t r 9r 
L x q: 

9Q 3F = 
9t 3x 


= S 


0 


( 12 ) 




The operators are applied in a symmetric way to avoid any numerical biasing of the solution: 

n + 2 ¥ t r r ..n 

q -L x L r L r L x q 

Each operator consists of a predictor and a corrector step. Each step uses one-sided differencing: 
Predictor: 


n+— 

q 2 =q"- 


At 

6Ax 


(7Fj - 8Fj_j + Fj_ 2 ) n 


n + 1 _ 1 


1 1 ^ 
n ^ ^ At / \ n H 

q +q 2+ ^( 7F i“ 8F i+i +F i^) 2 


Again, to avoid biasing, the sweep directions are reversed between operators: 

n n +2 T + -r +— T ~+L “ + a n 

q ~ L x L r r x q 


(13) 


(14) 


(15) 
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At the edges, flux quantities from outside the computational domain are required for the spatial differences. These 
quantities are computed using third order accurate extrapolations from the interior of the domain: 


F i + 1 =4F i -6F i _ 1+ 4F i _ 2 -F i _ 3 


( 17 ) 


MATHEMATICAL FORMULATION OF BOUNDARY CONDITIONS 

The primary objective in CAA is to accurately compute the oscillatory flow field. Since the computational domain is 
finite, numerical boundary conditions must be imposed at the grid boundaries. Improper specification of these boundary 
conditions results in artificial disturbances which can render the computed oscillatory flow solution unacceptable. This 
problem is most pronounced at subsonic outflows. As such, outflow boundary conditions for subsonic flows have received 
special attention and will be discussed herein. 


THOMPSON (1987) 

Thompson 6 ’ 7 describes how to decompose hyperbolic equations into wave modes of definite velocity and then how to 
specify boundary conditions for the incoming waves. The starting point of Thompson’s analysis is the nonlinear Euler 
equations. The essence of his approach is that one-dimensional characteristic analysis can be performed by considering the 
transverse terms as a constant source term. The 1-D characteristic analysis makes clear which waves are propagating into 
and out of the computational domain. The amplitude of the outward propagating waves are defined entirely from the 
variables inside the computational domain, while the amplitude of the inward propagating waves are specified as boundary 
conditions. For nonreflective boundary conditions, the amplitude of the inward propagating waves are set to zero. 


To illustrate Thompson’s approach, let us consider the nonconservative axisymmetric Euler Equations written in 
cylindrical coordinates: 


where 


Q t +AQ r +BQ x +C = 0 


(18) 



(19) 
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and 




C = 


1 

r 


-pv 

0 

0 

2 

-c pv 


(22) 


We first consider the radial boundary. The axial derivatives are grouped with the source term to rewrite Eq. (18) as: 

Q { + AQ r + K = 0 (23) 

Following Thompson’s derivation, the matrix A is diagonalized using similarity transformation. The eigen vectors are 
obtained by solving 


det(A — A.I) = 0, 

where I is the identity matrix. The left and right eigenvectors of A are calculated using 


(24) 


lJA = Vi T 
A^Vi 

l.T r . = 8-. 
i J U 


(25) 


The matrix S is formed, where the columns of S are the right eigenvectors r^ and the similarity transform is defined: 


A = S _1 AS = 


Xj 0 0 0 

0 0 0 

0 0 x 3 0 

0 0 0 X. 


(26) 


Applying this transform to Eq. (23) gives: 


S -1 Q t +AS -1 Q r +S -1 K = 0 


(27) 


or 


S -1 Q t + L + S -1 K = 0 


(28) 


This leads back to 


Q t +SL + K = 0 


(29) 


where: 


L. =A,.1. t Q 

i l l ^r 


( 30 ) 
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In this case, the eigenvalues are: 


= v-c 

X 2 = v 

= v 

X 4 = v + c 


( 31 ) 


and the left eigenvectors are: 



' 0 ' 


[ c 2 ] 


O' 


O' 

»!=• 

0 

> 1 — ^ 
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> ’ 1 3“' 

1 

, i — , 

0 

-pc 

’*2 

0 

0 

’ *4 

pc 


1 


-1 


0 


1 


Eigenvalues 1 and 4 are the velocities of sound waves moving in the negative and positive radial direction; eigenvalue 2 is 
the velocity for entropy advection; and eigenvalue 3 is the velocity at which vorticity is advected in the radial direction. 

Evaluating Eq. (30) gives 


L, =(v-c)(p r -pcv r ) 

L 2 =(v)(c 2 Pr - Pr ) 

L 3 =(v>K) 

L 4 =(v + c)(p r + pcv r ) 


Eq. (28) may be written as: 

Solving (34) for the time derivatives, we 


c 2 p t -p t - L 2 -c 2 K 1+ K 4 -R 2 
u t =-L 3 -K 2 = R 3 
P t + pcv ( = -L 4 - K 4 - pcK 3 = R 4 


(33) 


(34) 


obtain 



( 35 ) 


To apply the nonreflective Thompson boundary conditions, it is necessary to determine which waves (Lj-L 4 ) are 
incoming or outgoing. For example, if the boundary in question is the outer radial boundary, and the flow in the radial 
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direction is outgoing and subsonic, then the L l wave is incoming (velocity = v-c <0) and the waves are all outgoing. 
At this point, the amplitude of the incoming waves are set to zero, giving the boundary condition: 



(36) 


and the time derivatives for the boundary points are calculated using Eq. (35). 

The boundary conditions for the axial direction are obtained in a similar manner. After grouping the radial derivatives 
with the source term, Eq. (18) becomes 

Q t + BQ x +N = 0 (37) 

The eigenvalues and eigenvectors of B are computed, and a new similarity transform is defined: 


This time, the eigenvalues are: 


M = T -1 BT = 


0 

0 n 2 
0 0 
0 0 


0 0 ‘ 
0 0 
M-i 0 
0 H 4 


Mq = u-c 

^ 2 = u 

^ 3 = u 
M-4 = u + c 


and the left eigenvectors are: 



' 0 " 


fc 2 ] 


[0] 


o' 

m l = - 

-pc 

0 

' ,m 2 

0 

0 

' ,m 3 = 

0 

1 

%m 4 =< 

pc 

0 


1 


-1 


0 


1 


(38) 


(39) 


(40) 


In the axial direction, Eq. (33) become 
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( 41 ) 


Mj =(u-c)(p x -pcu x ) 
M 2 =(u)(c 2 p x -p x ) 

m 3 =<“(%) 

M 4 =(u + c)(p x +pcu x ) 


and Eq. (37) may be rewritten as: 


P t - P cu t = -Mj - N 4 + pcN 2 = Ej 

c 2 P t -P t =-M 2 -c 2 N 1+ N 4 =E 2 (42) 

v t =-M 3 -N 3 =E 3 
P t + pcu t = -M 4 - N 4 - pcN 2 = E 4 


The boundary conditions are evaluated as before, and the time derivatives are computing using: 



u 


t 



v 


t 



P,=}( E 1 +E 4) 


(43) 


GILES (1990) 

• Giles 8 derived boundary conditions based on Fourier analysis of the linearized Euler equations (LEE). Consider the 
unsteady part of the flow to be small disturbances, denoted by superimposed on a constant mean flow U, V. The 
axisymmetric LEE can be written as: 


where 


Q' t + A Q' r +B Q' x +c = 0 


(44) 


Q = 


(45) 



0 

V 
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0 0 


P 0 
0 0 

V I 
p 

c 2 p V 




(46) 
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and 


'U p 0 

0 U 0 

0 0 U 

0 c 2 p 0 


0 

J_ 

P 

0 

u 



-pv 

0 

' 0 ' 

-c 2 pv 


( 47 ) 


(48) 


An overbar denotes a mean flow quantity. 

Consider the disturbance to be in the wave form: 


Q'(x,r,t) = q R e i ( kx + lr - COt ) (49) 

where q R is a constant vector. Giles’ basic idea is that the dispersion relation for the above linear equation can be modified 
to prohibit propagation for waves with group velocities directed into the computational domain. This dispersion relation is 
nonlinear, which implies that the boundary condition must be a nonlocal one to satisfy this condition exactly. To produce 
local boundary conditions, the dispersion relation is expanded in a Taylor series around the one-dimensional solution. 
Boundary conditions for various degrees of approximation can thus be constructed. 

In the axial direction, the characteristic variables are defined as: 

Tj = p' -peu' 

Tj = p' -c 2 p' (50) 

r 3 = P' +p cu' 
r 4 = pev’ 


« 


* 


For a subsonic outflow condition (x=b), the r, characteristic is obtained using the Giles analysis; all others are ob- 
tained from the inner solution. The second-order Giles outflow condition is: 


3r, 


1 _ 


at 


= -u 


^4 + l4 
3r r 


-V- 


ari 

9r 


For a subsonic inflow condition (x=a), the fourth-order Giles boundary condition becomes: 



— V- 


dr 


2 3 _ 

dr at 




(c-U) 

2 


dr r 


-V- 


d r 3 

3r 


(c + U) 9r 3 (c - U) dt) 
2 3r 2 3r 


(51) 


( 52 ) 
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Likewise, for an inflow at an outer radial boundary (r=rmax), the characteristic variables are: 


c \ ~ P' -pcv' 

c 2 = p’ -C 2 p' 
c^ = p’ +pcv' 
c 4 = pcu' 


The subsonic inflow condition becomes: 



(c~V) 8c 4 

2 3x 3x 


9c 


3c„ 


— = — U- 
3t dx 


9 c 4 _ 9 c 4 (c + v)9c 3 (c-v)9c 1 

9t 3x 2 3x 2 3x 


( 53 ) 


(54) 


TAM AND WEBB (1993) 

Tam and Webb, 9 Bayliss and Turkel, 10 Hagstrom and Hariharan, 11 and Enquest and Majda 12 considered boundary 
conditions based on an asymptotic analysis of the Linearized Euler Equations with constant mean flow in the x-direction: 


“', +W x + f P', = ° 

v'.+UV +ip’ =0 (55) 

l X p I 

p', + Up' x+ pc 2 ( U ' x+ v ' r+ ^] = ° 

The above system of equations can be reduced to the convective wave equation for the pressure. For an outgoing 
wave solution at large distances, the boundary condition can be stated as: 


1 

V(0) 


P't+P' 


R 


P 

R 


(56) 


where 


R = Vx 2 +r 2 


V(0) = c 


X X, 

— M + 
R 



( 57 ) 


and 
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( 58 ) 





The above form of the pressure boundary condition was obtained by Tam and Webb, but one can show that it is 
identical to the pressure condition obtained by Bayliss and Turkel and others. 

However, Tam and Webb performed a rational asymptotic analysis for the full set of Euler equations. They showed 
that while the asymptotic pressure is only acoustic in nature, the velocity and density fluctuations contain both hydrody- 
namic and acoustic contributions. Using this, they arrived at a complete set of outflow boundary conditions which can be 
written in the form: 


= -Up' x+ ^(p- t+ Up' x ) 
c 

( 59 ) 

“', = - U "'x-=Px 

( 60 ) 

V t = - W x-=P'r 

( 61 ) 

= -V( 0 )[— p’ +— p' +— 

w Lr Kx r r rJ 

( 62 ) 


This boundary condition is applied at the outflow boundary; at the inflow boundaries, acoustic radiation boundary 
conditions may be used: 


<? t = -v(e) 


£<j +— q 1 + q ’ 
p v x R v r 


R 


( 63 ) 


where: 


Q' = 


( 64 ) 


TEST PROBLEM 

The test problem studied was that of a monopole in a uniform flow in the axial direction. The analytic solution is: 


p' (x,r,t) = 


J! [ 1+ ' I'm 2 -— ) 

4ttD D ) 


co cos(cdt) - sin(cot) 
D l 


( 65 ) 


u' (x,r,t) = 


4kD 


] / x ^ X 

— ~ — I M (ocos(cox) T sin(cox) 

P 2 cV D ) D z 


( 66 ) 


11 


Sr 


v'(x, r, t) = T 

4nD z 


cos(cox) - L- sin(cor) 


(67) 




( 68 ) 


where: 



P = a/i-m 2 

D = Jx 2 +(pr) 2 


T = t- 


D-Mx 

CP 2 


(69) 


For the described computations, 


n 

co = — 

2 

S = .01 (70) 

M = 0.5 


The computational grid starts at r = 0.5 and ends at r=24.5 in the radial direction. In the axial direction, the grid covers 
from -12 < x < 12. The grid has 300 x 300 points, and is not stretched in either direction. Tests were also performed with a 
finer grid of 450 x 450 points; while the overall errors decreased, the relative performances of the boundary conditions 
were the same. 

The point source is located at the origin. The analytic solution is specified as the initial conditions, and a CFL of 0.5 is 
used. With no mean flow, this gives 25 points per wavelength, and 191 time steps per cycle of oscillation. 

For all computations, the analytic solution was specified at the inner radial boundary. The upstream axial boundary 
and the outer radial boundary were both treated as inflow boundaries, while the downstream axial boundary was treated as 
an outflow boundary. 


IMPLEMENTATION OF BOUNDARY CONDITIONS 


Thompson 

The Thompson boundary condition was implemented in four ways: two were second-order accurate in time, while the 
remaining two were first order accurate in time. The axial outflow boundary will be used to illustrate the various imple- 
mentations of the Thompson boundary condition. 

Due to the split operators used by the code, the implementation of the Thompson boundary conditions was greatly 
simplified. The axial operator used by the code is defined by Eq. (12) as: 


3Q + 3F 
3t 3x 


= 0 


(71) 


% 


4 
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Thus, Eq. (37) becomes: 


Q t +BQ x =0 (72) 

since the radial derivatives and the source term do not appear in the axial operator. Using these definitions, Eq. (42) is 
rewritten as: 


p t ~ pcu t =-M i =u i 

c2 P t -P t =- M 2 =U 2 (73) 

v t =-M 3 =U 3 

p t +pcu t =-M 4 =U 4 


Using Eq. (73), the amplitude of the characteristic waves U r U 4 can be computed using either time or space 
derivatives. 

At the comers, the boundary conditions for both grid faces are applied. 

Thompson split/time (Th2t) . The first implementation of the Thompson boundary condition used the time derivatives 
to compute the amplitudes of the characteristic waves. To preserve the code’s second-order time accuracy, the boundary 
conditions were applied at the end of both the predictor and corrector sweeps in each direction. 

Thompson splii/space (Th2s'> . This variant used the spatial derivatives to compute the amplitude of the characteristic 
waves. The spatial derivatives were computed in the same manner as shown in Eq. (14) and (15), using extrapolation when 
required. Note that the derivatives are computed using the primitive variables, not the conserved quantities. Again, the 
boundary conditions are applied at the end of the predictor and corrector sweeps in each direction. 

Thompson/space CTh 1 s') . This version also used the spatial derivatives to compute the amplitude of the characteristic 
waves. This time, however, the boundary conditions were only applied at the end of the corrector sweep, using fourth-order 
accurate one-sided differences to evaluate the spatial derivatives at the old time level. First-order accurate time derivatives 
are used to update the flow variables at the boundaries. 

Thompson/time (Thlf) . The time derivatives were used to compute the amplitude of the characteristic waves. Again, 
the boundary conditions were only applied at the end of the corrector sweep, using first-order accurate time derivatives. 

Giles 


Two implementations of the Giles boundary conditions were performed, the first was first-order accurate in time, while 
the second was second-order accurate in time. 

Giles CGI) . In this version, the Giles boundary condition is implemented at the end of each time step, after both the 
radial and axial operators are applied. The spatial derivatives are computed using fourth-order accurate one-sided differ- 
ences, with the flow variables from the old time level. The time derivatives are computed using first-order backward 
differences. 

Giles split (G2i . Next, the Giles boundary condition was split into axial and radial operators. In this formulation, the 
boundary conditions are applied at the end of each predictor and corrector sweep, returning second order accuracy in time. 

The splitting procedure for the inflow boundary is described as an example of the methodology used. Eq. (50) gives 
the characteristic variables: 


13 


( 74 ) 


ij = p' — pcu' 

r 2 =p'-c 2 p' 

= p' +p cu' 
r 4 = Pcv’ 

The total Giles boundary condition for the axial inflow is given by Eq. (52): 






(c-U) ^4 , r 4 
2 ^ 3r r 


A 

/ 



-V 

V 



\ 

) 


(c + U)d r 3 (c-U) 3r ] 
2 3r 2 3r 


(75) 


Equation (75) is split into axial and radial operators, to be applied during the axial and radial sweeps, respectively. 
Since there are no axial derivatives, the boundary condition that is applied during the axial sweep is: 


3r, 


r =0 

3t 


dr. 


, 1 = 0 

3t 


3r 


^- = 0 

3t 


(76) 


During the radial sweeps, the boundary condition is: 


3r 3 

~3t" 


^i = _V^2 

3t dr 


(c-U) 


/3l 4 + i4 


dr 


^ 3r, 

dr 



^dr T ; 


(c + U) 3r 3 (c-U) 3r, 

2 dr 2 dr 


(77) 


The spatial derivatives were calculated from the primitive variables employing the differencing method used by the 
inner code. Extrapolation of the primitive variables was used when necessary. 

The axial outflow and outer radial inflow boundary conditions are formulated in a similar way. » 

Both Giles boundary condition formulations were unstable when applied to comer points. Since the performance of 
the boundary conditions at the comer points wasn’t a major area of interest, the analytic solution was specified at the comer # 

points for these calculations. To ensure that this was not causing errors of its own, a special comer treatment based on 
rotated 1-D characteristics was also tested; the performance of the Giles boundary condition was insensitive to the choice 
of comer treatments. 
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Tam and Webb 


In their original implementation of this boundary condition, Tam and Webb used a ‘boundary region’ of two points 
instead of one at each edge of the grid. No extrapolation was used in computing the spatial derivatives; this required con- 
structing stencils which only used variables at grid points in the computational domain. This is described fully in Ref. 9. 

To remain consistent, the Tam and Webb boundary condition was coded using only one boundary point. For comparison, a 
two point ‘boundary region’ variant was also tested; the results are equivalent. Two implementations of the Tam and Webb 
boundary conditions were performed. 

Tam and Webb (TW1 ) . In this version, the Tam and Webb boundary condition is implemented at the end of each time 
step, after both the radial and axial operators are applied. The spatial derivatives are computed using fourth-order accurate 
one-sided differences, with the flow variables from the old time level. The time derivatives are computed using first-order 
backward differences. 

Tam and Webb split (TW21 . Next, the Tam and Webb boundary condition was split into axial and radial operators. In 
this formulation, the boundary conditions are applied during each predictor and corrector sweep, retaining second order 
accuracy in time. The spatial derivatives were calculated from the primitive variables employing the differencing method 
used by the inner code. Extrapolation of the primitive variables was performed when necessary. 

Since the Tam and Webb outflow condition only corrects the pressure, the original code is used to compute the other 
flow variables. The pressure is computed using: 


p' t =-v(e) 


x , 

R P * 


(78) 


for the axial sweep, and: 


P't 



(79) 


for the radial sweep. 

For the inflow conditions at the other boundaries, the radiation boundary conditions are used. These are split as: 


Q’ 


t 




x 


for the axial sweep, and: 



(80) 


(81) 


for the radial sweep. 

Q’ is defined in Eq. (64). 


RESULTS 

Four points were chosen in the acoustic field to evaluate the various boundary conditions at. The first point (PI) was 
located at (x,r) = (-8.1, 20.5). The second point (P2) was at (8.1, 20.5). These points illustrate the effect of the boundary 
conditions on the comers of the flowfield, where the waves are not aligned with the boundaries. The third point (P3) was 
located at (0, 20.5). With no mean flow, this point shows the effectiveness of the radial boundary conditions when the 
wave velocity is normal to the boundary. The fourth point (P4) was located at (10.4, 2). With no mean flow, this point 
shows the effectiveness of the axial boundary conditions when the outgoing wave velocity is normal to the boundary. Also, 
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an axial line (LI) was defined at r = 12.5 to better illustrate the effect of implementation differences on the performance of 
the various boundary conditions. Figure 1 shows the location of these points and lines in the flowfield. 

At each point, the maximum amplitude of the pressure during each cycle was computed and compared to the 
analytic solution. The relative pressure error was used to compare the results. It is defined as: 


Again, the inner radial boundary for all computations was specified from the analytic solution. 

No mean flow fM = 0.0) 

Initially, the code was run for the case of a monopole in a still medium. All quantities are identical to the mean flow 
case with the exception of M = 0. While this is truly a radiation boundary, it was specified as an outflow boundary when 
using Thompson or Giles formulations. 

Figure 2 shows the maximum pressure contours for the exact solution. 

Thompson . The Thompson boundary condition was the easiest to implement, and was very robust. The only version 
that proved unstable was the implementation which used first-order time derivatives to compute the characteristics (Thlt). 
This version was always unstable due to comer problems, and was abandoned. 

All three implementations were basically equivalent. Figure 3 shows the maximum pressure contours for the Th4t 
implementation. Figure 4 shows the performance of the various implementations on the LI line. Figure 5 shows the 
performance of the Th4t implementation at the four test points. 

PI and P2 show the symmetry of the test case. The errors are the highest at this point, with values of nearly 25%. P3 
returns the lowest error, with 3% relative pressure error. At the lower comer, P4 shows 7% relative pressure error. 

The difference in the PI and P3 errors is due to the angle of incidence of the acoustic waves impinging on the bound- 
ary. When the wave velocity is normal to the boundary, the Thompson boundary condition performs best, as shown by the 
low errors at P3 and P4. 

Giles . The Giles boundary condition was the next boundary condition that was implemented. This condition was 
relatively easy to implement; it was also the most complex of the three due to its mixture of time and space derivatives of 
the primitive variables. The Giles boundary condition was stable everywhere except at the comers; since special treat- 
ments for the comer points were not given, the analytic solution was specified at the outer comer points. 

There was little difference between the performance of the G1 and G2 implementations of the Giles boundary condi- 
tion; Figure 6 shows the maximum pressure contours for the G1 implementation. Figure 7 shows their performance on the 
LI line. Figure 8 shows the results from the G1 version at the four test points. 

The Giles outflow boundary conditions also did not perform well for this test case, with relative pressure errors at PI 
and P2 of nearly 30%. At P3, the error is 1 1 %; while at P4, the error is less than 0.8%. The effect of the fully specified 
comer points on the far field error was not fully investigated. 

In other tests, a test case run in which the Giles inflow boundary condition is used instead of the outflow boundary 
condition gave much better results, with maximum errors of less than 7%. 

Tam and Webb . The Tam and Webb outflow boundary condition was the last formulation to be investigated. This 
condition does not have a special comer treatment; the comers are computed just as the other boundary points are. The 
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Tam and Webb boundary condition was stable everywhere, and relatively easy to code. Figure 9 shows the maximum 
pressure contours for the TW2 implementation. The second order implementation (TW2) was slightly better than the first 
order implementation (TW1). This is illustrated in Figure 10. Figure 11 shows the results of the TW2 run; errors at all 
four test points were less than 2%. 

Test runs were made using the radiation inflow boundary condition; performance was equivalent to that of the TW2 
outflow boundary condition. 

Mean Flow (M = 0.5) 

Figure 12 shows the maximum pressure contours for the exact solution. 

Thompson . All three implementations were all nearly equivalent. The two second-order versions (Th2s and Th2t) 
returned identical results, while the stable first-order implementation (This) was virtually indistinguishable. Figure 13 
shows the maximum pressure contours for the Th2t implementation. Figure 14 shows the results on the LI line from the 
three implementations. Figure 15 shows the results of the Th2t implementation at the four test points. 

The Thompson inflow condition did not perform well. For all three versions, relative pressure errors of 70% were 
computed at PI. It was theorized that the low performance of the inflow boundary condition in this case was due to the 
errors caused by the outgoing acoustic wave velocities not being normal to the boundary. 

The outflow condition performed much better; at P2, the relative pressure errors were below 20%. Here, the wave 
velocities are nearly normal to the boundary. 

The inflow condition specified at the outer radial boundary (P3) resulted in relative pressure errors of 20-30%. 

The Thompson boundary condition performed best at P4; the velocity of the outgoing waves were nearly normal to the 
boundary, and the proximity of the fully specified boundary kept numerical noise to a minimum. 

In general, the Thompson boundary condition was the least accurate of the boundary conditions tested. 

Giles . While the differences in the results were minor, the second-order Giles implementation (G2) slightly outper- 
formed the first-order (Gl) version. Figure 16 shows the LI results from the two implementations. Figure 17 shows the 
results of the G2 run at the four test points. Figure 18 shows the maximum pressure contours for the G2 run. 

The Giles inflow condition worked very well, with 

relative pressure errors below 3% at PI and P3. At P4, the proximity of the fully specified boundary kept the relative 
pressure errors below 5%. 

The Giles outflow condition did not perform as well as the inflow condition, with the error at P2 nearly reaching 20%. 

The Giles boundary condition is the second best of the three formulations tested. 

Tam and Webb . Figure 19 shows the results of the runs using the Tam and Webb boundary conditions; again, the 
second order implementation (TW2) was slightly better than the first order implementation (TW1). Figure 20 shows the 
results of the TW2 run at the four test points. Figure 21 shows the maximum pressure contours for the TW2 run. 

Radiation inflow conditions were used at the upstream axial and outer radial boundaries; these resulted in errors of 5% 
at PI, and 3% at P3. 

The Tam and Webb outflow condition performed very well, giving errors of 1% at P2 and P4. 

The Tam and Webb outflow boundary is the best of the formulations tested, resulting in lower errors than the exact 
boundary conditions. The only drawback to this formulation is its present restriction to uniform mean flows. 
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CONCLUSIONS 


Several boundary conditions for computational aeroacoustics were evaluated, namely, Giles (1990), Tam and Webb 
(1993), and Thompson (1987). The same discretization scheme was used to compute sound radiation from a monopole 
placed in a stagnant or moving stream, but with the boundary treatment as described by one of the aforementioned 
approaches. For each boundary condition, various implementations were tested to study the sensitivity of their perfor- 
mance to the implementation procedure. In general, slight improvement in the performance of each scheme can be 
achieved via the implementation procedure. 

With the best achieved performance for each boundary treatment, the following conclusion can be drawn regarding 
their comparative performance: 

(A) For outflow boundary treatment, the only acceptable performance was that of Tam and Webb. The performance of 
the other schemes might be acceptable only in special cases wherein the flow is nearly one-dimensional, perpendicular to 
the boundary. But even under such conditions, their performance did not surpass that of Tam and Webb. 

(B) For inflow boundary treatment, Giles boundary condition was acceptable, as well as the radiation boundary 
condition of Tam and Webb applied to inflow treatment. The Thompson inflow boundary treatment resulted in consider- 
able reflection near the inflow boundary. 
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Figure 1 . — Computational domain. 
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Figure 2. — Maximun pressure contours for uniform flow about a monopole 
t (M = 0.0; analytical solution). 
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Figure 3. — Maximum pressure contours for uniform flow about a monopole 
(M = 0.0; Thompson 2nd order/time). 
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Figure 5— Uniform flow about a monopole (M = 0.0; Thompson 2nd order/time). 



Figure 6. — Maximum pressure contours for uniform flow about a monopole 
(M = 0.0; Giles 1 st order). 
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Figure 7. — Uniform flow about a monopole (M = 0.0; Giles). 
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Figure 9. — Maximum pressure contours for uniform flow about a monopole 
(M = 0.0; Tam and Webb 2nd order). 
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Figure 10. — Uniform flow about a monopole (M = 0.0; Tam and Webb). 
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Relative pressure error 




Figure 1 2. — Maximum pressure contours for uniform flow about a monopole 
(M = 0.5; analytic solution). 
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Figure 13. — Maximum pressure contours for uniform flow about a monopole 
(M = 0.5; Thompson 2nd order/time). 
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Figure 14.— Uniform flow about a monopole (M = 0.5; Thompson). 
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Figure 16. — Maximum pressure contours for uniform flow about a monopole 
(M = 0.5; Giles 2nd order). 
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Figure 17. — Uniform flow about a monopole (M = 0.5; Giles). 
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Figure 1 8. — Uniform flow about a monopole (M = 0.5; Giles 2nd order). 
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Figure 19. — Maximum pressure contours for uniform flow about a monople 
(M = 0.5; Tam and Webb 2nd order). 
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